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Abstract 

Spectra of circularly polarized harmonics is calculated by numerically solving the Time- 
Dependent Schrodinger Equation for a 2D model of Ne atom using circularly polarized funda¬ 
mental with counter-rotating second harmonic laser fields. We demonstrate strong asymmetry 
between left- and right- circularly polarized harmonics when a ground state with p-type symmetry 
is used. It arises due to the circular polarization of individual attosecond pulses in the generated 
pulse train. Reducing the length of the counter-rotating drivers and introducing a small time-shift 
between them allows to generate a single elliptically polarized attosecond pulse. 
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High Harmonic Generation (HHG) in atoms and molecules is a highly nonlinear process 
which up-converts intense infrared laser field into the extreme ultraviolet (XUV) and soft 
X-ray radiation W- The emitted light can be used to track quantum dynamics underlying 
the nonlinear response or as a table-top source of bright, coherent, ultrashort pulses 

[M2]. 

In the latter case, generation of circular or highly elliptic high harmonics and/or attosec- 
ond XUV pulses is very important. Such pulses would find numerous applications, e.g. in 
chiral-sensitive light-matter interactions such as chiral recognition via photoelectron circular 
dichroism (PECD) (see e.g. [T31 - fT5] ). study of ultrafast chiral-specific dynamics in molecules 
(e.g. P| Hj), and X-ray Magnetic Circular Dichroism (XMCD) spectroscopy (e.g. [IEH21]), 
including time-resolved imaging of magnetic structures (e.g. [IBH22])- Table-top sources of 
sub-100 fs, or even attosecond, chiral pulses would be a real breakthrough for laboratory- 
scale ultrafast studies. Not surprisingly, search for schemes enabling the generation of short, 
coherent XUV pulses with tunable polarization is a very active area of research, see e.g. 

PIESH3H]. 

Importantly, the control over polarization is desired not only for individual harmonics, 
where it has just been demonstrated [35, dB], but also for individual attosecond pulses, 
both isolated and in a train, where robust and practical scheme is still lacking. We show a 
way to solve this problem, proposing a practical scheme for the generation of highly elliptic 
attosecond pulses, both single and in a train. 

An elegant solution to generating individual high harmonics with circular polarization 
has been found by W. Becker and coworkers [28, 29| 31, 32j. It relies on combining circularly 
polarized fundamental field with a counter-rotating second harmonic. The resulting electric 
field peaks three times within one cycle of the fundamental, producing three ionization 
bursts. The electron promoted to the continuum near the peak of the instantaneous field can 
successfully revisit the parent ion within about half-cycle, emitting an attosecond radiation 
burst [281 [29]. 

This approach has now been very successfully used in [35, 36j, demonstrating generation 
of bright, phase matched high harmonic radiation. Importantly, tuning the ellipticity of 
one of the fields allows to tune the ellipticity of the generated high harmonics from linear to 
circular [35 j. While the theoretical interpretation of this control is a matter of debate [331 (35] , 
the approach is very promising. However, until now the possibility of extending this scheme 


2 


from controlling the polarization of individual harmonics to controlling the polarization of 
isolated attosecond pulses looked far from straightforward. 

Indeed, the driving field dictates that the direction of electron return rotates by 120° 
three times per cycle. Consequently, recombination with an s-state yields three linearly 
polarized attosecond bursts per cycle, with polarization rotating by 120° from burst to 
burst [29j. Thus, while each harmonic is circularly polarized, the same does not apply to 
their superposition. 

This can also be seen in the frequency domain. The harmonic lines are at energies 
{n-\-l)uj-\-2nuj = (3n + l)o; and na; + (n+l)2a; = (3n + 2)a\ In centrally symmetric medium, 
and for circularly polarized driving fields, the selection rules dictate that the f2 = (3n + l)u 
line has the same circularity as the fundamental while the fl = (3 n + 2 )cu line has the same 
circularity as the second harmonic, (12 = nuj + 2 nut = 3 nuj is parity forbidden) [HU E2I ESI 
3(T . 40i [41]. Thus, the harmonics have alternating helicity. Adding harmonics of alternating 
helicity with equal intensity yields an attosecond pulse train where each subsequent pulse 
has linear polarization rotated by 120°, in concert with the time domain picture. 

Suppressing every second allowed harmonic line, e.g. Q (3n+2) = (3n + 2)cu, would solve 
the problem of generating individual attosecond pulses with circular polarization. O. Kfir 
et. al. [3f)J suggested that such suppression can be achieved by optimizing the phase¬ 
matching conditions in gas-filled hollow fiber and reported substantial suppression of the 
lines 12( 3n+2 ) = (3n + 2)ui. 

First, we show that relative intensities of the counter-rotating harmonic lines strongly 
depend on the orbital momentum of the initial state. For an initial p-state (as for a Neon 
, Argon, or Krypton gas), the harmonics co-rotating with the fundamental field can be 
much stronger then those co-rotating with the second harmonic. The effect is found with 
the contribution of both degenerate sub-levels, p+ and p—, included in the calculation. As 
a result, circularly polarized attosecond pulses are generated already at the microscopic, 
single-atom level, see Fig. [l] Additional help from phase-matching is a bonus, but not 
necessary. 

Second, we extend the scheme to generation of isolated attosecond pulses. We show that 
when the counter-rotating driving pulses become relatively short, e.g. 7-8 fs for the 800 
nm driver and its second harmonic, one can generate an isolated attosecond pulse, or a 
controllable train with 2 or 3 pulses, by tuning the time delay between the fundamental and 
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the second harmonic. 

To demonstrate these effects, we numerically solve the time dependent Schroedinger equa¬ 
tion (TDSE) for a 2D Neon-like model atom, for counter-clockwise (+) polarized fundamen¬ 
tal and clockwise (-) polarized second harmonic. We show that the harmonics generated 
from orbitals with m= ±1 differ from those generated from s orbitals in two important 
ways. Firstly, the height of the adjacent left- and right- circularly polarized harmonics can 
differ by an order of magnitude, with m= 1 state favouring harmonics co-rotating with fun¬ 
damental and m= —1 state favouring harmonics co-rotating with 2u> field. Secondly, once 
the two contributions are added coherently, + polarization continues to dominate in a broad 
spectral range, leading to highly elliptic circularly polarized attosecond pulse train already 
at the single-atom level. Our findings are in accord with [36] (see Ne spectra in Fig. 3 of 
[36] ). where such disparity was attributed to phase matching. 

We solve the (TDSE) in the length gauge (atomic units are used throughout unless stated 
otherwise): 


d . . 


T + V[r) + r • E(£) <b(t, r) 


(1) 


The 2D model potential is taken from [32] 


Z(r) 

where Z(r) = 1 + 9exp(—r 2 ) and a = 2.88172 to obtain the ionization potential of Ne atom 
I p = 0.793 a.u. for the 2p orbitals. The Is state has an energy E\ s = —2.952 a.u. and the 
2s energy is E 2s = —0.217 a.u.. For reference calculations we use Is as the initial state but 
keep the same ionization potential taking Z(r) = 1 and a = 0.1195. 

The laser electric field is 


E(t) =E ir ■ f(t) ■ (cos[off] + cos[2a;t])a: 

+E ir ■ f(t ) • (sin[off] — sin[2uff])y (3) 

where f(t) is the trapezoidal envelope with 2 cycle rising and falling edges and 5 cycle 
plateau (in units of fundamental). The cn-field rotates counter-clockwise (+). The second 
harmonic rotates clockwise (-). 

The TDSE is propagated on a 2D Cartesian grid using Taylor-series propagator with 
expansion up to 8th order [43] . A complex absorbing potential 

V c (x) = rj ■ (x - x 0 ) n (4) 
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laser frequency 

u) 

0.05 

(A = 911 nm) 

laser electric field 

E 

±^ir 

0.05 

(1= 0.88 • 10 14 W/cm 2 ) 

grid step size 

dr 

0.2 


time step size 

dt 

0.005 


propagation time 

T 

1250 

(30.2 fs) 

maximal grid extent X max 

±60 


absorbing boundary xq 

±36 



TABLE I. Parameters of the calculations in atomic units unless stated otherwise. 

with i) = 5x 10 -4 and n = 3 is used to avoid non-physical reflections from the boundary. 
Other simulation parameters are summarized in Table [T| 

Convergence was tested with respect to the absorbing potential, the time step and the 
spatial grid. Note that HHG in bicircular fields is dominated by very short trajectories [2BJ. 

The initial wavefunctions were obtained using imaginary time propagation filtering out 
the ground state wavefunction to obtain p x and p y orbitals. The p± states are defined as 
p± =Px ±i p y The laser intensity was kept such as not to exceed 5% ionization and to 
avoid strong shifts and mixing of the degenerate atomic orbitals described in [42]. The 
spectra were obtained by performing the Fourier transform of the time-dependent dipole 
acceleration, evaluated at every 0.5 a.u. 

The results are robust with the variation of the pulse length, the shape and length of 
its rising and falling edges, laser intensity, and wavelength: we performed calculations from 
A = 600 nm up to A = 1200 nm. 

Fig 1(a) shows reference spectra obtained for Is initial state of the model potential with 
I p of Neon. It agrees well with previously published results [2.8, [321 ESI EE], the harmonics 
come in pairs ( n + 1)oj + n2u> = (3 n + l)m and nuj + (n + l)2u; = (3 n + 2 )u> of similar heights. 
The left harmonic in the pair has the same polarization as the fundamental field, the right 
harmonic follows the 2a j driver. The harmonics 3nu> are parity forbidden. 

Figures [D 3 and [lj; show spectra for the p + and p_ initial states. For the p + initial 
state, the harmonics that have the same polarization as the driving IR field are preferred. 
For the p_ initial state, the harmonics with the same polarization as the 2u> driver are 
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stronger. There are additional spectral variations in the plateau region, different for p + and 
p_ orbitals. There is also a qualitative difference between the below-threshold (<I p ) and 
above-threshold (>I p ) harmonics, showing that the evolution of the photoelectron in the 
continuum is critical for the observed propensity in the harmonic strengths. 

Figure [l]d shows the spectra obtained from adding the contributions from the p + and p_ 
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FIG. 1. Spectra for (a) Is, (b) 2p + , (c) 2p_ initial states and (d) equal mixture of 2p + and 2p_ 
states. Colors mark harmonics co-rotating (red) and counter-rotating (blue) with the uo field. 
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orbitals coherently, as required. In the plateau region, harmonics with the same polarization 
as the driving IR field dominate over those with opposite polarization. 



Time (fs) 


FIG. 2. Gabor Transformation spectrogram of the harmonic spectra intensity and ellipticity from 
Is (a) and (b) and 2p (c) and (d) orbitals. Colour in (b) and (d) indicates the ellipticity of 
the spectral components in the regions of the spectrogram where the amplitude of the spectra is 
significant.The horizontal dashed lines mark the l p . 


The sub-cycle dynamics of the emission process was analyzed using the Gabor Transform 
(GT) [41J of the time-dependent acceleration dipoles a(t): 

GT[Q, t 0 } = — [ dto(t)e-^e- (t - to)2/(2T2) (5) 

2n J 

where we have chosen T = 1 /3cn. The reference spectrograms for the Is initial state in 
Figures^ and [2]) show the time-dependent intensity (a) and ellipticity (b) for time-resolved 
spectra, in the regions where spectral amplitudes are significant. As expected, there are 3 
radiation bursts per oj cycle with linear polarization, as predicted in [28, 29, [32] . 

Figures [2[: and HP show the same spectrogram for the 2p state, i.e. the coherent su¬ 
perposition of the radiation from p + and p_ states. Although the signal strength in the 
spectrogram is similar to the s orbital, the ellipticity of the emitted radiation is very differ¬ 
ent. Three distinct regions can be identified: (i) below threshold region, where the ellipticity 
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is mostly negative; (ii) the middle region, where the ellipticity is high and positive and (iii) 
near cutoff region where the emitted radiation is mostly linear. The energy region (ii) of the 
spectrogram coincides with the spectral window in figure [ljl where the difference between 
clockwise and counter-clockwise harmonics is the greatest. 



Time (fs) 


FIG. 3. Time-resolved XUV emission from a 2p orbital, for time-delayed 8 fs base-to-base 800 
nm and 400 nm pulses, (a) Spectral intensity and (b) time-dependent ellipticity for the perfect 
overlap of the two-pulses, (c) Spectral intensity and (d) time-dependent ellipticity for the two-pulse 
delay of 2 fs shows that single, highly elliptic attosecond pulse is generated. Colour indicates the 
ellipticity of the spectral components in the regions of the spectrogram where the amplitude of the 
spectra is significant.The horizontal dashed lines mark the l p . 

Application of bi-circular fields is naturally extended from the generation of an attosecond 
pulse train to the generation of an isolated attosecond pulse, using short driving pulses and 
changing the time-delay between them. Indeed, the harmonic emission driven by circular 
fields is only possible when the two counter-rotating circular pulses overlap. Given high 
nonlinearity of the overall process (including ionization), it will be limited to the temporal 
window where the two driving pulses overlap with nearly equal and high intensity. This idea 
is tested in Fig. [3| which shows time-resolved spectrograms and ellipticity of the emitted 
light for A = 800 nm and 400 nm counter-rotating drivers with base-to-base duration of 8 fs 
















and sin 2 envelope (3 u field cycles), for two different time delays. In case of perfect overlap, 
3 attosecond pulses are generated. Delaying the low-frequency driving pulse by 2 fs (| of u 
field period) yields a single attosecond pulse with highly elliptic polarization. 

What is the physical origin of the HHG sensitivity to the angular momentum of the 
initial state? The energy and angular momenta that the electron accumulates from the 
laser field while propagating in the continuum is transferred to the harmonic photon upon 
recombination. The matrix elements associated with recombination are the complex con¬ 
jugate of the photoionization matrix elements. In 2D one photon ionization with the field 
co-rotating with the initial state is much more likely than with counter-rotating field. This 
is a direct analogue of Fano-Bethe propensity rules fff5j and is also the case for Rydberg 
states co-rotating and counter-rotating with the field [4H, 1 7]. 

Consider the harmonic spectra from p + orbital. The right circularly polarized harmonics 
result from the (n + 1)uj + n'2uj pathway. The recombination step is conjugated to photo¬ 
ionization from p + state with a co-rotating field, favoured by the propensity rules. The 
left-circularly polarized harmonics result from the noj + (n + 1)2uj pathway. The recombi¬ 
nation step is conjugated to photo-ionization from p + state with a counter-rotating field, 
dis-favoured by the propensity rules. This explains the relative heights of the harmonic pairs 
for the p + initial state. The same analysis explains why harmonics co-rotating with 2u> field 
are preferred for the p- initial state. 

But why is p + dominant over pJ? The answer lies in the stronger effect of the lower- 
frequency (counter-clockwise) field on the continuum electron, which leads to higher popu¬ 
lation of the continuum states with positive angular momentum than with the negative one. 
The more probable recombination from such states is to the p + state, by emitting light with 
counter-clockwise polarization. 

Finally, we comment on the importance of the carrier-envelope phase (CEP) stabilization 
in this scheme. As long as the relative phase between the two pulses, u> and 2u>, is locked, 
changing the CEP will rotate the polarization ellipse of the attosecond pulse but will not 
alter its high ellipticity. This property, in combination with the possibility of using relatively 
routine durations of the two driving pulses, makes the scheme extremely attractive for 
practical implementation. 
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